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‘The ofreutt shown in Pie, 1 above, is a ‘lew pase filter suitable 
fer suppresatac random noise appearine in the Inout wolbage a(t), GCon- | 
sider the tape voltage e(t) as esmpesed of two parte: 
| ‘Inch dostrabte signal input veltace of constant level (0.¢.)* 
ie noise, whieh varies rather rapidly with tine, 
fy a suitable chetes of resistance and capacitance 14 1a possible 
or reduce the noise voltage sopesring in the output v(t) to « negligible 
Shed, whtle the desired #ignal is transnitted intact. 
Movewor, tise varying efgnals present in the input signal in adéstion 
te the constant devel: wilh — filtering to @ smaller or greater 


Tt ony perhaps seca rather vii te vowed a D.C, voltage or voltages 
varying Linearly with tine as signal voltages; one is accustaned to 


: in ne se meen oe However, = 


pe antes re ree ee ee | 
| y tn antonatte: Af 


extent alongside with the noise voltages, depending upon the rapidity of 
time variations 

in order to obtain an optimus conpromise between the two conflicting 
requirements of mexiaum ecoureey in signal reproduction and minimum noise 
output, it is sometimes necessary to vary the cireuit paranetere, The 
analysis of the clreuit, after the variable parameter has envuumed its 
final values can eraelly be accomplished by conventional cireuit theory; 
however, fow investigations have been made in the past concerning the 
behavior of output voltege during transition time, in fret, past analyses 
hove often borrowed the solutions valid fer constent parameters, replecing 
the conatant perameter by ite timeevarying equivalent. Thie approximate 
method of selution very often yields surprisingly scourate reaulte; 
however, it breaks down in other oases, 

it is the aim of this paper to investigete the exact performance of 
the reeistanceecapacitance series circuit with time worying resistance, 
and to cheek the range of error involved in the approximate meted of 
evalueting the cireuit behavior. The resistance variation to be con= 
eidered primarily is a linear variations higher order terme are considered 
in the appendixes 


Figs 2 


The differential equation governing the above circuit is 

(1) OR(t) dy(t) + v(t) = o(t) 
oh linear sakes of resistance is exsumed; Le@e 

(2) A(%) s Ro + Ryt 
For spavenianes define 

(5) OR(t) = OR, + GAYE ET, it 

Two types of input voltages e(t) are to be studied, namely power 

‘ae containing finite numbers of terms, and sinusoids, In many 
applications 1t is sufficiently accurate te regard the power series ae 
morely a straight line or parabolic variation; emphasie will be pleced 
on these two special cases, However, extension of the solution to any 
| desired finite number of terme ie given in the appendix. 
| Yince the resistance variation dincusced in this paper is 2 tine 
function only, equation (1) is a linear equation; consequently the supere 
. ponition principle may be applied. Thus, one may consicer each term of 
7 a(t) by itself separatelyy one solve for MPe, terms of v(t) whieh are due — 
to this individual term of ‘o(t) tones Ae a final step the solutions 


obtained in thie manner may be added, resulting in the solution due te all 
of the terms of e(t) combined. 

Linearity and superposition are very convenient in thie study, since 
they permit separate investigation of power series and sinusoidal signels. 
Furthermore, one ic assured thet neglection of higher order terms in e(t) 
will not affect the solutions of v(t) corresponding to the lower order 


terms of e(t). 


ate 


The pertinent form of equation (1) is 


t @, for constant signel voltage 
(4) an ” i) oe, + v(t) = i + © i* for pe god tignal voltage 


e, + eit +e, for parabolic signal voltage 

Solutions to these differential equations are presented in Table i, 
subject to the boundary condition that et t ge © the condenser is charged 
to the conetant input voltage e,3 at t = 0 simltaneous variations in 
resistance and input signal begin. This boundary condition approxim tes 
the conditions of the practical applications of this probleme 

The solutions for in 2 O are also presented, since the previous 
epproximete solution, presented in Table Ii, has merely replaced - by 
t+ Tt wherever |, appears in Table I, Uols 2, 

The signel error €(t), defined as v(t)-e(t) ie tabuleted in Table ILI 
ap obtained by both the exact and approximate methods. 

E(t) corresponding to e(t) = e, + e,t is of immedicte interest in 
design problems; inspection of Table Lil reveals thet the signal error is 


independent of e 5 furthermore for the case of e(t) = e, + ,t the signal 


l 
error is directly proportional to @° 

For plotting purposes it will be most convenient to define a new 
variable 5 =.= &(t). $ may perhaps be termed “per-unit signal error". 


Sinee the poe of E(t) and 6, are volts and volte per second reste 


i 
pectively, the proper dinedsienk of $s are seconds. 

Pes caloulated by both the exact end approximste hethodi 18’ ptéttes 
against time with te and Te as parameters. The following ranges of 
variables and parameters are presented ae most nearly representative of 


automatic fire control applicationa: 


l--Time: 0-20 seconde 


Za 15 seconis 

ae =) to 3 | 

: The fully dram ourves with index Bex! nepreeadt solutions 
‘Aine by the exeet inethed, while broken Line curves with index 

Ze ; e ap" "saponins solution is obte ‘ined: by the sppreximate method, A 
"singe curve with index. * "Pex ard, Zam" indicstes that the differences 


i mn between the tito ji eiealins ordinates as ent culated by the two 


: methods are less than .1 nit of B(.1 second); this fieure become 


ne _ tomparable to the finite curve thickness; consequently 3 plot of the 


a approximate ourve was not warranted, PRON SHARES 3h become 


- mecessary to obtain more exact figures, especially in the neighborhood 


of & ® § seconds, where differences. between Rapp and Bex become 
, t ‘ “ ‘ Ti 


Re “Anherently very email, one may resort to the curve entitled " Sapp < 


: Rex VS. qi ah & = 5 seconds", i caimdcimeats there are caccified in 


; the, ratio ef 101. Se ea 


| The pereenbaze error in Renee voltage oe St as tabulated in 
‘Tables IV, ¥. for r elt) | = &, +e t at Ae 10 “seconds end 20 seconds ves- 


ial pectively. 


oe ot <8 eng ea ie it is ‘povaible 
is et aa dal ~ t A oe ay a tae 
Nie * 53 pene ae =: 


: te ‘simi taneous ‘tabulate | Sow’ one eevee and one negative ratio 


: of at constant te The oe (=) ie in the tabulations refers to 
a es outare HESS “we ‘the “ of the sie pe vice versa. 


Table I- Exact Solutions 
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Table TT- Signal Error £(4 
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Table IV 
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Per cent Error in Output Voltage Neca at te lO sec. 
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Table IV- Continued 
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Table V 
Per cent Error_i utput Voltage | Ot t=O sec 
Die core to CRM Sh ee 


| 716.361 +4.04| =10.20 
Poiehs PROD TS 12 1240 x (26), 7) 100 
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Oe clusions 


in wash gn work the following consideratione may be ant are 


l-~Absolute magnitude oi error £ (t) = v(t) ” e(t) 


2--Percent magnitude of error d s neat 
P 7 e s 


fwo methods of calculating €(t) and Tare available, namely the 
exact and approximate methods, 

For linear signel voltage and linear resistence varistion, the 
graphical and tabuler results offer quick means to compute 2(+) and Hy . 

To find€(t) corresponding to given l, ’ is end time, it is merely 
necessary to mltiply the appropriate ordinate by “e)¢ This may be done 


graphicelly by adjusting the scale of ordinates to the required @)3 Te- 


calling that $2 =&t) « 
e 
The absolute sdpistate of the error becomes greater with time, !, 


ry and 3... For positive [| the approximate absolute error becomes 


i. 
progressively greater than the exact error, ae each ef the aforementioned 
paremetere become largers however, the opposite is true for negative r. 

Thus the maximum error corresponds to Te = 5 seconds, Pe = 35 and 
+ = 20 seconds for which &ex(20) = 8.2 & volte, Capp(20) * -10,5 @, volts, 
a difference of -2.1 ®, volts. On the other hand a more accurate plot of 
thecdifference at t = 5 sec, shows that Capp(5) = 2ex(5) is a3 8) volts 
corresponding to i = 5 secs, by = 3. 

The percent error "| depends on the additional parameter eg » In- 

| 

epection of Tables IV and V indicates that for a wide range of 2 at 
~~ 


~q- 
& oh fr. 
> ery oe 
4 


Hypa 
aH cee 
% a eae * 


+= 10 secs t = 20sec, The difference between “fapp and 7 ex does not 
exceed 5, To reduce || to & winiwum @ chould be highs 


€ 
Sc eg L : 
Ranges of linear signal voltages not plotted or tabulated ae well 


as saduheiks voltages may be investigated with the aid of the expressions 
of Tables I, II, ill. It i9 hoped thet eufficient preliminary investigation 
hos been made concerning higher order voltage and resistence terms, as 
presented in the Appendix, The following limitations arise in the con- 
| siderations of higher order terms: : 
i la-The mumber of parameters increases to e poini where general ealeulee 
tions ean very cumbersomes ror any given set of values, however, 
~ ealeulations are feirly simple. 

2-siane of the ieiwidens obtained there are of necessity presented in. 
integral forms it would probably be necessary to perform these integrations 


on the differential analyzere 


i , =. 
-? 


7 


ie mana ae at ide 


4g) 


Xe * + fe erp 


‘es aah mer, Ue ayiliemy ated 


=e tebch te ieee soni: ets oat 


oy ACR 


ut ra HM torenrtal acetbed fans nine clas ache 
Soe Roger ad went Meme 


zt 
wt 4B vinwk) ir, +74)? 


"hi 


1 
(t 
mye 


7 ) ep wh 
eo) mh CEse HH 


4 wt ae.) one ws eek aca Te LT 4" MS ite Pp ope, : 
Sawer ote Letham! @yeeitiag on one Gelieriag fexe, t ; - 
tT a { : : 
! {tec | rs 
4 . : 4 1; 2. ‘ 7 . ‘. ~~ A i? 
gi » ‘ert, Ore we, mw dae wt. ve 38 H of. 4, 
{4 } ‘ ee F (es) ‘> ey we free eee OT) Sie aR ce ecm latin ri i i AG oat et hy 
i f “yy act * 
: } by Ths 
; uw, 9 ed 
| 7 : 7 
‘~ 
‘irae 8 8th = ‘ Meee %¥ 4m rn ty r atta re er % m dupe my 7 
4 Dit Cel OORT oe ‘ woe Ee OGL Woy id a Des Nb 9 Mo Prd {oe see 
a 
i “1 ah Bethea sae ye ie aeturaie si = ‘Lae Wea as Pepe * i, ms a ite aah _ : 
“igh ait ‘fh Yul a oars yi or ia “4. i He 2) Trad sf ae < ‘ Hey pre ah aS ot Gt i ey pO 6ORD. ing” 
3 
he ad 
P af ay fae o 8 “5 - bc . * i Moses ' ae 4” % 
4, ¥ Wigs FAY ad WF Ps Py) 1S 3, L Ad ee fe PM Ge ty.) a2 yre! ae oid Ty - ohn a 2 “') 
ie " > : 
‘04 n , we Geo det ea eae ‘ oe ee ae he we A eee ge, * [ee 4 
ope * Be mPa kl ye wd mae UA Beni | a8 ‘ # % he Ws 2s aa tin: dhepyne wut 
‘a 
eee ui " 1% yr ts J 7 wy & = * ee 3 "y * 
rari aoe, Fh Soe oxo heralater? Fog |i BB a 
% atin Oe ee m eae fee 
ee Te ut ; ? om Chee ake 
' 
| “i ithe ‘ 4 v2 re De paMeoge, 
ug he a Gee 4 as whew Twin le 
“her 
Pde Ne Th MAR eS we CO Ba a Opera 
ih 
‘De vey Ak 
are a 
Bits wa Os uty nit, be i? 
om PRL Va oe F 


“li 
Sglution ef the Differential Equation 


—_ wate): of the fom T, + T, t, and excluding the case for 
tT, = 0 which has been solved elsewhere in this paper, the solution 
of the differential equation for sinusoidal signal volteces requires 
the solution of the following integral equations 


f . _2 a ‘ 
(5) Vv +: ~ fe eos wt + B sin ot) (% ryt)? al (t,, $hyt) Ts. 


We ean set the arbitrary constant K, equal to V(o) in which case we 
have our integral equation in the following form: 


: ; Tv 
3 Tv a t a 
(6) V=| Vo) + / ‘ obese, 8 8 fort, i 
log -) 


1, +t, ¢ 


When the exponent ey ~ 1 is a small positive integral number, we are 
able to evaluate _ integral directly by the process of expansion of 
the term (t ottyt) tT - eae Selutions are tabulated for Tt *1, 2, 35 
and 4, in addition to the case in which Tt, = © which is the constant 
parameter case, The solutions are tabulated for the general case in 
_ Which A and B have non-zero vernent this is the case for any phase 
angle of the composite input weve, Solutions are also tabulated an 


S the case in whieh A 20, corresponding to a phase angle of — and 


‘for the ease in which B= 0, corresponding to a phase angle of zero. 
For higher peeitans integer values of be » for negative integer 
values of ; and for all fractional values ot recourse must be 


had to a differential analyzer, a Fourier coefficient analyzer, 

or numerical hand computation, In attacxing the problem, however, 
we are particularly interested in the range -.3<7y <.3 and it is 
to be noted that the tabulated solutions for 7) 0 (special case) 
and 7, = .25 fall in this range. Further, the case wherein 


us = .33 lies at the extreme positive end of the range. 


oo] Jee 
é ‘ Differential ‘ equation: ; ave) ae W(t) = o(t) 
PR ayaa he at GA(t) ~~ Gre) 


Signal voltage:  e(t) = A cos ot +B sin ot 


| Case I: CR(t) = constant # Te 
CR(t): : 
Case IIs R(t) = i tt 


Steady state solution 
(7) V steady pie ire wi AaB 
for Case I: State = smi COSUt + sinwt 
w@ 22 +1 wT 2 +1 


In the solutions which follow for Case II, the solutions are 
tabulated for the general case, for the case in which A = 0, and 
for the case in which B= 0, 


Since e(t) = A coswt + B sinwt: 


e(t) = a2 4 pe cos (w t-6) 


where @ = taa ~" 


Hence, in the case where A = 0, 


6s ¥ and e(t) = B sinut 


In the case when B = 0, 


@ = 0 and e(t) =A cosut 


For © any value the general solution applies, in response to the 
signal voltage e(t) = A cos ot ¢ B sinwt, | 

In order to indicate the variation of V as X function of ™ 
for small positive salves of 7y» the accompanying curves have 
been plotted under the following conditions: 


= 0, 1/3, 1/2) 3 
AmBei 

es 

fo #1 


The solutions then become? 


20: Vesint 

4 21/31 ¥> a [: + (t*99t413) sin t = (t24ute2) cos :| 
fi aif ve ane [2 + (3e2t) sin t = (14) cos ‘| 

A 21: vs te [2 tein t = one | 


For Fi equal to sere, which is the case when CK(t) = constant 
= ts we havg a constant coefficient differential equation to solve, 
the signal veltace being a sum of sinusoids which may resolved into 
a single sirmsoid with an initial phase angle. As is to be expected 
the solution for the response is e constant amplitude sinusoid, As 
a 1 takes on small positive integer values the solutions becone more 
complicated but, as ean be seen from the equations and the curves, 
the response ecain involves a sinusoidal form, though now we have 
a vertical displacement of the sinusoid due to the aprearance of 
the transient term which was neglected in the ease where a = 0, 

In this latter case it was considered that at the starting of 

problem time the transient had already died out and steady state 
conditions existed, For small positive integer values of % the 
sinusoidal term ere modified by amplitudes that are uot constant 


whe 


... that are funetiions of time. The net result is that the 
| = "responses to sinunedéal. inpate, in the case of CR(t) variable, 
: are of 2 damped sinusoidal form with the damping factor being 
- variable and not following a definite pattern. In general, 
ree Hs - inereasing A decreases the amplitude of the response (en 
senepiden is seen in the ease of 7 = 4 where the initial 
swing of the response is to a value greater than unity) though 
te make this eouparigon it is necessary te take into account. 
the phase shift ie iepanstine of the resultant cycle introduced 
“by inereasing ae from. its nero value. Pearwasing 7% towards 
sero has the effect of changing the doletion towards the constant 
| moplitude agmanpca form, : 


i Hy el give ‘an Andication of the effect. ati changing frequency 
is the solutions have been plotted when only one input sinusoid 
is present for two valnas of angular frequency, =] and w= 2, 
est other parameters have the same values for the two solutions. 
The curves show that as the frequency is inereased the irregularly 
| damped sinusoidal fom responses are dncsemnaid in amplitude and 
die out sooner, 
Again confining our attention to sinusoidal signal voltages 
‘the question comes up as to whether or not we are justified in 
substituting the setual variation of CR(t) inte the solution for 
cCR(t) = constant, and to whet extent error is introduced by this 
preeedite, } | | 
Our Anput signal voltage is e(t) = A coswt + B sinvt, 
To obtain an average indication, consider a given solution in 


=] be 


two ways; first, let A = 0, B =1 which eorresponds to 6 = - 
and second, let A = 1, B = 0 which corresponds to 6 = 0. Denoting 
the first solution by Vp and the second by V,, we will plot 

vA" * Vp" for both the exact method of solution and for the 
approximate method wherein we substitute the actual variation of 
CR(t) in place of CR(t+) = constant, Considering all other 


parameters fixed and the same for all solutions, we have: 


For TY alex 

Approximate t u,4 r ¥,* - + 
(let), 1 

penis %," 7 V,° “ 1.5 + sin t = cos ¢ 
(1 ++)? 

For 7 =: 

Approximate: V 2 Vv," = saeonsinstidcieincen 

Pi A B (Lo 5 t)2 43 


(1 +.5t)? 4 5 + gin t ~ (1 +.5t) cos t 


Exact * 4“ ¥," = 
(1 + .5t)4 


The curves show that the responses die cut more slowly as 
Y, is decreased toward zero, and that agreement between the exact 
and approximate solutions is better for smaller 7. That the 
response shovld die out more slowly as % is redueed is to be 
expected since when q = 0, the sum of the squares of V, and V, 
will be equal to a constant. The disparity between the approximate 
and the eect solutions as a function of T} is another question. 
Obviously, at some value of 7 and for lower values, we will be 
within the limits of error specified by a given design problem 
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“ Clasotest nathods were preferred to the nethod of Laplace 
| raneforontions in the solutions to the differential equations, 
" siinee the transforne corresponding to the varying circuit 
_ paraneters are net generally avai lable in the literature, 

Ae an exanle the foveal solution for V(+) eorrespending 
- te-a parwbolte olga) wi.) be presented here! 

: pn (4) should perhaps more accurately road: 
2 ae ae emt) _ 
- Because. of the Linearity of this equation one mar salve 

SOUS sug ita ame betig 0 and ¢, setec, # 

ix soameay. barat tue haa of the solution are conmonly tered 


* Vit 
(t) = ry Re 2 +0 


ren tay, : steady state  sattons, 
a "Por the transients 
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; oe where K is an arbitrary conotant . 
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Por the steady state solution 
Assume that 

% ) ¥ = t 2 
(8b) (+) A + 3% + Vt 


4,2. V(t) is assumed to be axpressible in a power series terminating 


in the same decree as ft) 


Differentiatine (%>) 
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(@e) at 1 * al ot 
Multiplying (8c) by (7% #%t) and adding to (8) there results 


Gt) (Xp +Ht) (Vy + B¥ot) Vy + vy + V7 


BR + | + Res + 2% | t + V0. 42%) ¢? 
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Therefore: 
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Solving (9b) simultaneously 
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eas 1H) (1s2h ) | boundary 
as condition. 


Solving (94) for ® 
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The method of "assuming" the steady state solution may be employed 
to obtain solutions to any number of finite terms appearing in the expansion 
for e (t), 1.e. let 


(10) e (t) Se,+eqttest* --~--- - e,t” 
Assume as a solution to equation (1) for CR(t) = rid. lg" 

T = 2 n 
(11) T(t) = VgtVit+Vot* +------- Vnt 


By the previous method one obtains as solutions 
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This scheme of "assuming" the steady state solution will succeed as 
long as the resistance is constant or varies linearly with time; however, as 
soon as parabolic or higher order terms are admitted into the expansion of 
R(t) the method fails. The reason for failure is easily evident: teal 


is an expansion whose highest degree term is of one degree. lower than the 
highest degree term of V (t). ial tiphying Ai by a term containing t* 
or t raised to a power higher than 2 will cause the degree of the left hand 


side of equation (1) to exceed the degree of the right hand side, making it 
impossible to equate coefficients. 


Oa Ob a, 7 


However, since equation (1) is a first order linear differential 
equation its solution can always be placed in integral form if not closed 


form. For our equation: 


(13) on (+) MeLivey =e (t) 


aV(t) , V(t) = Sih 
(U4) a ath GR(t 
For the homogeneous equation (right hand side = 0), we have: 


Vit o 
as a + Ud, 20 


Separating the variables we obtain: 


at = ee a 
CRC) t 


[> fi 
log -* - hiss 


Se 
~ CR(t 


(16) = 7) he 
When the right hand side of the equation is not equal te zero, we assume the 
same form of solution but in this case the constant x is replaced by u, 


where u is now a time function. 


Then: fates 


Y= we at : dt 
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t dt Gait) 
! a ay © 
ay = du fa 
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Substituting into the original equation we have: 
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Then: 
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én M t . OR(t 
dt CRI 
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fen(e) 
Therefore: 
is 
~ | Gare 
V= we - 
(19) tee - [Secs 
CR (t CR(t 
v= Ko + sl. dt e 


This method of attack, involving the solution of an integral equation, is 
used in this paper when (1) R(t) is assumed to consist of terms giving 
rise to a variation greater than linear, and (2) the signal voltages are of 
sinusoidal form. 

For R(t) consisting of terms of higher order than linear subse- 
quent analysis shows that the "transient" part of the solution can be placed 
inte closed form, but the steady state solution cannot. However, an integral 
equation in itself constitutes appreciable progress compared to the differ- 
ential equation form, 

R(t) = Rot R, t+ Ryt* 


(19) | Consider 
CR(t) = To4 Tt 4Tot® = T2 (t®+ + t+ a2) 
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In so far as sinusoidal signal voltages are concerned, to 
| illustrate the method of obtaining the results tabulated previously in the 
paper, the following sample calculations are shown for the case in which 


Tote 
ck = z 71 
=. + se e(t 
ty CR(t) 


(21) Let e(t) = A cos u t4#B sin u t 


cR(t) =I 4+ Ty t 
avy v - A cos ut+B sin ut 
a igi 7, Te ety * 
Then: ; at a dt 
v= ro fsoe ns Bsinut ,/ Toy * wl) oJ fot tit 
— Spe tree eee ee cr ee Ca 
Tot Ti * | 
A cos ut i. fa log ("} 0+ Tt) - log (T-+Lt) 
ot /y 
+-1 _ 
corm 
(22) re figt co von a neiye 7) a (To47,%) : 
Since T, = 1: 


[ ae Ts 
r=/ mt] a cos u t +B sin u t) af Pe ius eae 


: A B ot 
(23) relate cin ut Zoos ut| ree 
To evaluate the arbitrary constant K, we must go back to the steady state solution 
for R constant and substitute t = 0, giving us V(0), our boundary condition 
for the variable resistance solution above. 


For constant resistance we have: 


Wy V = Acosutts sinut 


- We 


The steady state solution is found as follows: 


Let Vo = steady state solution for V 


Let Vo = M cos uttN sin ut 


a e- uM sin uttu N cos ut 
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Squating sine terms and cosine terms separately: 
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(25) V2- = cos u t+— 2 sin u t 
a aed ul +h 
For the variable resistance case, V(0) = Vo (0). 


(26) Volz A- uTy B 
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(27) K3 = lg MEN ~ 
“ara: 


The complete solution for V in the variable resistance case is then: 


= (A-uT, B) 
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